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ABSTRACT 

A dormant supermassive black hole lurking in the center of a galaxy will be revealed when a star 
passes close enough to be torn apart by tidal forces, and a flare of electromagnetic radiation is emitted 
when the bound fraction of the stellar debris falls back onto the black hole and is accreted. Although 
the tidal disruption of a star is a rare event in a galaxy, « 10~^ yr~^, observational candidates 
have emerged in all-sky X-ray and deep ultraviolet (UV) surveys in the form of luminous UV/X- 
ray flares from otherwise quiescent galaxies. Here we present the third candidate tidal disruption 
event discovered in the Galaxy Evolution Explorer (GALEX) Deep Imaging Survey: a 1.6 x 10*'^ 
erg s~^ UV/optical flare from a star-forming galaxy at z = 0.1855. The UV/optical spectral energy 
distribution (SED) during the peak of the flare measured by GALEX and Palomar Large Field Camera 
imaging can be modeled as a single temperature blackbody with Tbb — 1.7 x 10^ K and a bolometric 
luminosity of 3 x lO'*^ erg s~^, assuming an internal extinction with E{B — V)gas — 0.3. The Chandra 
upper limit on the X-ray luminosity during the peak of the flare, Lx(2 — 10 keV)< 10*^^ erg s~^, is 
2 orders of magnitude fainter than expected from the ratios of UV to X-ray flux density observed in 
active galaxies. We compare the light curves and broadband properties of all three tidal disruption 
candidates discovered by GALEX, and flnd that (1) the light curves are well fitted by the power- 
law decline expected for the fallback of debris from a tidally disrupted solar-type star, and (2) the 
UV/optical SEDs can be attributed to thermal emission from an envelope of debris located at roughly 
10 times the tidal disruption radius of a « IO^Mq central black hole. We use the observed peak 
absolute optical magnitudes of the flares (—17.5 > Mg > —18.9) to predict the detection capabilities 
of upcoming optical synoptic surveys. 

Subject headings: galaxies: nuclei — ultraviolet: ISM — X-rays: galaxies - black hole physics 



1. INTRODUCTION 

Stellar dynamical models predict that once every 10'^ — 
10^ yr a star in a galaxy will pass within the tidal dis- 
ruption radius of the central black hole, Rp < Rt — 
i?*(A/BH A^*)^^^, and will be torn apart by tidal forces 
(Magorrian & Tremaine 1999; Wang & Merritt 2004). 
After the star is disrupted, at least half of the debris is 
ejected from the system, and some fraction (0.1 — 0.5)M^,, 
remains bound to the black hole and is accreted (Rees 
1988; Ayal et al. 2000). The fallback of debris onto the 
black hole, and stream-stream collisions between pre- 
cessing debris orbits (Rees 1988; Caimizzo et al. 1990; 
Kim et al. 1999) produce a luminous electromagnetic 
flare that is expected to peak in the ultraviolet (UV) /X- 
rays (Ulmer 1999). 
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Detections of tidal disruption events are of interest 
because they are a signpost for the presence of a dor- 
mant supermassive black hole in the nucleus of a galaxy 
that is otherwise quiescent and undetectable (Frank & 
Rees 1976; Lidskii & Ozernoi 1979; Rees 1988). The ob- 
served properties of tidal disruption flares have the po- 
tential to probe the mass and spin of supermassive black 
holes in distant galaxies, independently of locally estab- 
lished scaling relations between the mass of the central 
black hole and its global host galaxy properties (Geb- 
hardt et al. 2000; Ferrarese & Merritt 2000; Graham et 
al. 2001; Marconi & Hunt 2003). Tidal disruption events 
may also reveal the presence of intermediate mass black 
holes in globular clusters and dwarf galaxies (Ramirez- 
Ruiz & Rosswog 2008). However, for the most massive 
black holes, the tidal disruption radius is smaller than 
the event horizon and stars are swallowed whole without 
disruption (for a solar- type star Merit ~ IO^Mq, Hills 
1975). Future wide-field time domain surveys that will 
coincide with the LISA mission have the exciting poten- 
tial to use the electromagnetic detection of a tidal dis- 
ruption event to localize the burst gravitational waves 
that is produced when the star approaches the tidal dis- 
ruption radius (Kobayashi et al. 2004), and to look for 
tidal disruption events from recoiling supermassive black 
holes after a binary merger (Komossa & Merritt 2008). 

The most convincing candidates for a tidal disruption 
event are luminous, UV/X-ray flares from galaxies which 
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have no ongoing active galactic nucleus (AGN) activity 
for which a flare could be attributed to the AGN's central 
engine. Strong candidates have now emerged in all-sky 
X-ray surveys {ROSAT All-Sky Survey (RASS), Don- 
ley et al. 2002; XMM-Newton Slew Survey, Esquej et 
al. 2006) and in the deep UV Galaxy Evolution Explorer 
(GALEX) Deep Imaging Survey (DIS; Gezari et al. 2006, 
2008). In a systematic search of 2.9 deg^ of GALEX DIS 
FUV and NUV data discovered two luminous UV flares 
from non active, early-type galaxies (Dl-9 aX z — 0.326 
and D3-13 at z — 0.3698) that were also detected in the 
soft X-rays in follow-up Chandra imaging. Because the 
DIS fields overlap with the optical CFHT Legacy Sur- 
vey (CFHTLS) Deep Survey, we were able to extract 
simultaneous optical difference imaging light curves of 
the flares in the g, r, and i bands with an excellent ca- 
dence of days. The broadband spectral energy distribu- 
tions (SEDs) of the flares, and their light curves were 
in excellent agreement with the theoretical expectations 
for emission from a tidally disrupted main-sequence star 
onto a central black hole of few times lO^Af© (Gezari et 
al. 2008). 

The new candidate presented in this paper was discov- 
ered in five years of GALEX FUV and NUV DIS obser- 
vations of the 1.2 deg2 DEEP2 23h field. A large am- 
plitude. Am ^ 2 mag, FUV flare was detected in 2007 
from a source that was not variable in the three yearly 
epochs prior. The source, hereafter referred to as D23II- 
1, was followed up with Keck LRIS optical spectroscopy, 
which identified the host as a star-forming galaxy. We 
identified D23II-1 as a candidate tidal disruption event, 
and followed it up with optical and X-ray imaging obser- 
vations nearly simultaneously with the peak of the flare, 
and months later, and with a late-time optical spectrum 
one year later. The paper is organized as follows, in Sec- 
tion [2] we describe the optical, UV, and X-ray imaging 
and optical spectroscopy of the flare and its host galaxy, 
in Section [3] we classify the host galaxy as a star- forming 
galaxy, and compare the SED and light curve of the flare 
with what is expected for a tidal disruption event, in Sec- 
tion |4] we compare the light curves and broadband prop- 
erties of all three GALEX candidates and compare them 
to all of the candidates reported from UV and X-ray sur- 
veys to date, and in Section [5] we conclude by using the 
observed optical properties of the flares to predict the de- 
tection capabilities of the next generation of optical syn- 
optic surveys. Throughout this paper, calculations are 
made using = 70 km s~^ Mpc~^, = 0.3, Oa=0.7 
and a luminosity distance for D23II-1 of = 900 Mpc. 
UT dates are used throughout, and magnitudes reported 
are in the AB system and corrected for a Galactic ex- 
tinction of E{B - V) = 0.035. 

2. OBSERVATIONS 
2.1. GALEX UV Photometry 

The GALEX FUV (Ae// = 1539A) and NUV (Ae// = 
2316A) observations of the DEEP23H_02 DIS field from 
2004 to 2008 were coadded into nine epochs with expo- 
sure times ranging from 6 to 17 ks. We performed aper- 
ture photometry with the center fixed to the centroid 
of the source in the epoch with the highest signal-to- 
noise detection, a =23h 31m 59^53, S = +00°17'14'.'57 
(J2000.0). We used a 6'.'0 and 3'.'8 radius aperture in the 



FUV and NUV, respectively, with aperture corrections 
from Morrissey et al. (2007) of 0.87 and 0.58, respec- 
tively. We chose a smaller aperture in the NUV in order 
to avoid contamination from a bright star located 45'.'8 
to the west of the source, and we added 0.19 mag of 
systematic error in quadrature to the photometric error 
to account for variable contamination in the aperture by 
the bright star. The sky background is assumed to be 
the pipeline-generated SExtractor catalog sky value at 
the location of the source (Berlin & Arnouts 1996) in 
a reference image constructed from a very deep coadd 
{tcxp — 100 ks). This avoids an overestimation of the 
sky background that has been found to occur in shal- 
low observations in comparison to deep coadded images 
(Morrissey et al. 2007). We determine the photometric 
errors by measuring the dispersion of the differences of 
~ 5000 sources in each epoch from the "true" magni- 
tude in the reference image as a function of magnitude. 
We also determine the mean offset between each epoch 
and the reference image by measuring the mean of the 
Gaussian distribution of the difference magnitudes for 
all sources brighter than 21.5 mag, resulting in offsets 
for each epoch from the reference image in the range of 
-0.039 to -f 0.058 mag in the FUV -0.023 to -1-0.020 mag 
in the NUV. Figure [T] shows the FUV and NUV GALEX 
images during the large amplitude flare, and Table[l]lists 
the photometry for each epoch. 

Figuredlshows the FUV and NUV light curve of D23H- 
1 and its FUV— NUV color over time. The flux of the 
source is constant within la in the yearly epochs from 
2004 — 2006, and the maximum peak-to-peak amplitude 
of the flare from its steady state on 2006 September 10 to 
its peak on 2007 September 29 is an increase in flux by a 
factor of 9.0 ± 1.5 in the FUV and 4.7 ± 1.2 in the NUV. 
The FUV— NUV color is 0.7 mag bluer during the peak of 
the flare, indicating that the SED of the source during its 
"flaring state" is different than during its "steady state" . 
The GALEX observations constrain the peak of the flare, 
to, to be tQ = 2007.74 ±0.07. 

2.2. Optical Photometry 

The optical counterpart to the flaring UV source in the 
DEEP2 DRl survey catalog from 1999 to 2000 (Coil et 
al. 2004) is a resolved galaxy with a half-light radius in 
the R band of 0'.'918 detected with B = 20.665 ± 0.005, 
R = 19.296 ± 0.001, and / = 18.825 ± 0.001, and in the 
Sloan Digital Sky Survey (SDSS) DR6 catalog from 2001 
September with u = 21.52 ± 0.25, g = 20.172 ± 0.027, r = 
19.281±0.018, i = 18.859±0.019, and z = 18.677±0.062 
(Adelman-McCarthy et al. 2008). We obtained optical 
imaging in the g band with the Palomar Large Field 
Camera (LEG) on UT 2007 October 15 (with 1'.'5 see- 
ing) and November 14 (with 1'.'9 seeing) corresponding 
to 16 and 46 days after the peak of the flare (or to + 16 
days and to + 46 days), with an exposure time of 1620 
s and 1800 s, respectively. The spatial resolution of the 
ground-based imaging is insufficient for determining the 
morphology of the galaxy beyond the fact that it is re- 
solved. Using the SDSS g-band images as a template 
we performed digital image subtraction using the High 
Order Transformation of PSF And Template Subtraction 
(HOTPANTS^), which is based on the original algorithm 

^ http:/ /www. astro. Washington. cdu/users/beckcr/hotpants. html 
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of lAlardI (|2000f) . In both epochs we find a residual point 
source consistent with the location of the host galaxy 
D23H-1 (Figure [3]). Using several nearby SDSS stars for 
reference, we measure g-band magnitudes of 22.25 ±0.12 
(2007 October 15) and 22.37 ±0.17 (2007 November 14) 
for the flaring source. 

2.3. Spectroscopic Observations 

We obtained an optical spectrum with Keck I LRIS 
and the 400/8500 grating on 2007 September 16 with a 1" 
slit resulting in a spectral resolution of AA — 5.7 A, and a 
wavelength coverage of 5700 - 8300A at to - 13 days, and 
with Keck II DEIMOS and the 1200 line/mm grating on 
2007 October 7 with a 1" slit resulting in AA = I.SA and 
a wavelength coverage of 4600 - 7300A, at io + 8 days, 
both with texp = 1800 s. The spectra are shown in Figure 
4. The LRIS spectrum was reduced with standard IRAF 
routines, and the DEIMOS spectrum was reduced using 
the DEEP2 DEIMOS data pipeline, which is optimized 
for data taken with the 1200 line/mm grating as used 
in this work. We calibrated the flux scale of the LRIS 
spectrum by integrating under the SDSS z-band filter 
curve and scaling the flux to match the i magnitude. We 
measure the redshift of the host galaxy from the narrow 
emission lines to be z = 0.1855. 

The narrow emission-line ratios and their la errors de- 
termined from Gaussian fits (shown in Figure [5]) to the 
spectra after subtracting a template for the stellar ab- 
sorption lines, are: [O III] A5007/H/3 = 0.6 ± 0.1, [N II] 
A6583/Ha = 0.41 ± 0.03, ([S II] A6716 + A6731)/Ha = 
0.22 ± 0.06, and an upper hmit of [O I] A6300/Ha < 
0.035, using a 3tT upper limit on the [O I] fiux measured 
from propogating the standard deviation of the contin- 
uum over 5 pixels. Figure [S] shows the line ratios in the 
standard diagnostic diagrams (Baldwin, Phillips, & Ter- 
levich 1981; Veilleux & Osterbrock 1987, Ho et al. 1997) 
and in comparison to the empirical dividing line be- 
tween star-forming galaxies and Type 2 AGNs measured 
from SDSS (Kauffmann et al. 2003) and the theoreti- 
cal boundaries between star-forming galaxies, Seyferts 
and low-ionization emission-line regions (LINERs) from 
Kewley et al. (2001). The line ratios fall within the 
boundaries for HII regions of star-forming galaxies in all 
three diagrams. In the [O III]/H/3 versus [N II] /Ha di- 
agram, however, the object lies very close to the bound- 
ary between star-forming galaxies and composite galax- 
ies from Kewley et al. (2001), and the boundary be- 
tween star-forming galaxies and AGNs from Kauffmann 
et al. (2003). This could allow for some contamination 
to a pure star- forming spectrum by an AGN nucleus. To 
investigate this, we plot a "mixing line" of a pure star- 
forming spectrum with a coordinate of (+0.45, +0.37) 
with a LINER nucleus with a coordinate of (0.0, +0.3) 
with gray dots in Figure [6l The position of D23H-1 is 
consistent with a mixture in which 90% of the Balmer 
emission is from a star-forming galaxy, and 10% of the 
Balmer emission is from a LINER nucleus. A mixture 
with a Seyfert spectrum at (+0.1, +0.6) would allow for 
only 5% of the Balmer emission to be from a Seyfert 
nucleus. In both cases, the LINER or Seyfert nucleus 
would contribute ^ 30% to the total [O HI] line emis- 
sion. Thus, although there may be some contribution 
from an AGN, the narrow emission lines are dominated 



by star formation. 

The spatial profile of the spectrum at the location of 
Ha and at the continuum near a rest wavelength of 6500 
A (which is dominated by host galaxy starlight) both 
have full width at half-maximums (FWHMs) that are 
equivalent to the seeing during the observation (esti- 
mated to be 0'.'8 - 0'.'9 from the FWHM in the spatial 
direction of the standard stars observed before and after 
the target). Thus, we do not have the spatial resolution 
necessary to determine if the Ha line is extended and fol- 
lowing the spatial distribution of the stellar continuum 
light (as would be the case if it was tracing star forma- 
tion) , or if it is more compact and coming only from the 
nucleus (as would be the case if it was powered by an 
AGN). 

Motivated by the fact that the narrow-line ratios clas- 
sify the galaxy spectrum as a star-forming galaxy, we fit 
the UV/optical SED measured by GALEX, SDSS, and 
DEEP2 in the steady state with a subset of 25 Bruzual & 
Chariot (2003 ; BC03) galaxy templates with solar metal- 
licity and a single burst of star formation with ages from 
3.2 X 10^ yr to 2.0 x 10^ yr. For these template fits, we 
have assumed that the UV and optical flux in the steady 
state is dominated by stars in the galaxy. We simultane- 
ously fit for internal extinction by adding extinction to 
the galaxy templates using the stellar extinction law from 
Calzetti (2001) with E{B - V)stars = OAAE{B - V)gas, 
for E{B - V)gas — — 0.6. We find the least-squares fit 
for templates which have a specific star-formation rate 
(SFR) of 1.68 X 10-i°Mq yr"! (M^/Mq)"!, where Mg 
is the mass of the galaxy, an age of 2 Gyr, and internal ex- 
tinctions of 0.17 < E{B - V)gas < 0.30. Figure H shows 
the spectra and the best-fitting BC03 galaxy template 
with E{B — V)gas — 0.3. The stellar absorption features 
in the spectrum, including Ca II, G-band, Mg lb, and 
Na ID, match the galaxy template very well, except for 
the Balmer series lines which are filled in from the con- 
tribution of narrow-line emission. Because the DEIMOS 
spectrum was not corrected for the response of the in- 
strument, we divide out the shape of the spectrum, and 
multiply it by the shape of the best-fitting BC03 galaxy 
template. 

Due to the possible contribution of a LINER nucleus 
at the 10% level or a Seyfert nucleus at the 5% level to 
the Balmer emission, it is possible that some fraction of 
the steady-state UV flux is non-stellar. If we repeat the 
same fits to the steady-state optical SED, and exclude 
the GALEX UV photometry, we can fit for lower extinc- 
tions, < E{B - V)gas < 0.30, lower specific SFRs of 
0.70 - 1.68 X 10-i°(Mg/MQ)-i, and older ages of 2 - 3 
Gyr. The best-fitting template to the optical data not 
allowing for internal extinction is the same as the best- 
fitting template to the UV and optical data allowing for 
internal extinction. Thus, the UV photometry is entirely 
consistent with the age and SFR fitted to the optical 
photometry as long as there is some internal extinction 
added. Therefore, instead of an excess of UV flux above 
the star-forming galaxy template that would suggest a 
nonstellar origin, some internal extinction in the galaxy 
is required to explain the observed deficit in UV flux. 

Note that there is an age-metallicity degeneracy, such 
that higher metallicity galaxies can produce redder op- 
tical colors that mimic an older stellar population. This 
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degeneracy can be broken by comparing the strength 
of the metal absorption-Une features, which increase in 
strength with higher metalhcity (see Figure 10 in BC03). 
Given the good match of the solar metallicity BC03 
template to the stellar absorption features in the spec- 
trum of D23H-1, invoking an older/lower metallicity or 
younger/higher metallicity template is not necessary. 

2.4. Internal Extinction 

The measured Ha/H/5 ratio is sensitive to the galaxy 
template chosen to subract the stellar continuum. If we 
use the best-fitted BC03 template to the broadband pho- 
tometry, and add the range of best-fitted internal ex- 
tinctions, 0.17 < E{B - V)gas < 0.30, this results in 
Ha/up = 4.1 — 4.4. These values are larger than ex- 
pected for case B recombination ([Ha/H/3]o = 2.87). 
Because the [O IIIJ/H/3 versus [N II] /Ha line ratio for 
D23H-1 is close to the boundary between star-forming 
galaxies and AGNs shown in Figure [SI it may be that 
D23H-1 has a "composite spectrum" , with some contam- 
ination by an AGN. In AGN spectra there is a contribu- 
tion to the narrow Ha line from coUisional excitation 
that increases the Balmer line ratio to (Ha/H/?)o = 3.1 
(Osterbrock 1989). The observed Balmer decrement is 
= (Ha/H/3)o/(Ha/H/?)ohs = 1.4 — 1.5 when assum- 
ing (Ha/H/3)o — 2.87, and 1.3 — 1.4 when assuming 
(Ha/H/3)o = 3.1, and thus (including all the systematic 
uncertainties) imphes E{B — V)gas = 0.3 ± 0.1 using Ta- 
ble 2 in Calzetti (2001). It is encouraging that both the 
BC03 template fits to the broadband photometry and 
the Balmer decrement indicate E{B — V)gas ~ 0.3. 

The observed Galactic extinction-corrected Ha lumi- 
nosity, L(Ha) = 8.3 X 10^° erg s""'^, translates to an un- 
extincted line luminosity of Lhq — (1.2 ±0.1) x 10'*^ erg 
s~^, and assuming that the Ha line is powered by hot 
stars, an SFR= (0.95 ± O.O8)M0 yr'^ (Kennicutt 1998). 
The B-band absolute magnitude of the galaxy is Mb = 
B - DM - K - Ab ^ -18.4 ± 0.2, where we have used 
the conversion B = Bab + 0.163, K = —1.24 measured 
from the galaxy template, DM — 51og(dL/10pc), and an 
extinction of Ab = 2.22E{B - V)gas = 0.7 ± 0.2 using 
the extinction curve for starlight from Calzetti (2001). 
This absolute magnitude corresponds to a galaxy mass 
from the best-fitting BC03 template of (4± 1) x lO^M©, 
and an SFR = (0.7 ± 0.2)Mq yr"\ which is in agree- 
ment with the SFR estimated from the unextinguished 
Ha luminosity within the inherent uncertainties in such 
techniques (e.g., at the factor of 2 level). 

2.5. Late-time Spectrum 

We obtained a follow-up spectrum 1 yr after the peak 
of the flare to look for high ionization He IIA4686 or 
iron lines that could appear as a light echo to a flare 
in the EUV/soft X-ray photoionizing continuum, as was 
reported to be the cause for fading high ionization lines 
in a broad-line AGN by Komossa et al. (2008). D23H- 
1 was observed with the Low-Resolution Spectrograph 
(LRS) and the Gl grism on the Hobby-Eberly Telescope 
(HET) on three separate occasions: 2008 November 3, 
17, and 18. On each observation two 1200 s exposures 
were taken through a 2" slit, oriented at the parallac- 
tic angle, with a final spectrum with AA = 14 A, and 
a wavelength range of 4090-10,830 A. The reductions 
were carried out in a standard manner. Correction of 



telluric absorption bands and (relative) flux calibration 
was carried out using the spectra of spectrophotometric 
standard stars observed on the same night and reduced 
in the same manner as the spectra of the source. 

The spectrum shows neither the appearance of new 
emission lines, nor any change in line ratios since the 
spectra taken in 2007. Since the spectrum is still dom- 
inated by the host galaxy hght, this suggests that (1) 
there is not enough gas within a light-travel time de- 
lay of 1 yr from the central black hole to reverberate 
the flaring continuum or (2) that the bulk of the flare 
emission did not extend into the EUV/X-ray wavelength 
range (hundreds of eV) that is necessary to excite these 
high-ionization lines. 

2.6. Chandra X-ray Observations 

We obtained 30 ks target of opportunity (TOO) ob- 
servations with Chandra ACIS-S on 2007 October 1-3 
{to + 3 days) and 2008 January 23 {to -1-116 days). Us- 
ing a circular aperture with a radius of 4 pixels and a 
background annulus with radii of 8 and 12 pixels, the 
galaxy has a < 2ct detection above the background from 
0.2 - 10 keV of (1.3 ± 0.8) x lO""* cts s-\ or a 3cr upper 
limit of < 2.4 x lO^"' cts s^^. This corresponds to an un- 
absorbcd flux density at 1 keV of < 1.7 x lO^'^^ erg cm~^ 
s^^ Hz~^, for an AGN power law with F = 1.7, and a 
Galactic absorption column density of Nh — 2.03 x 10^° 
cm~^. The total unabsorbed flux in the 0.2 — 2 keV 
band is < 8.5 x 10^^^ erg cm~^ s~^ which corresponds 
to L(0.2 - 2keV) < 7.2 x 10^^° erg s-\ and in the 
2 — 10 keV band is < 1.1 x 10~^^ which corresponds 
to L(2 - 10keV)< 9.3 x 10^° erg s'^ On 2008 January 
23 the galaxy was not detected above the la level with a 
3(7 upper limit of 2.1 x 10^** cts s^^, corresponding to an 
upper limit to the flux density at 1 keV of < 1.4 x 10"'^^ 
erg cm~^ Hz~^, and an upper limit to the 0.2 — lOkeV 
luminosity of < 1.4 x 10*^ erg s~^. If we assume that the 
X-ray emission suffers the same amount of extinction as 
the gas in the HII regions, then for Nh = (1.7±0.6) x 10^^ 
cm~^, the X-ray luminosities get multiplied by a factor 
of 1.3-1.5. 

3. INTERPRETATION 

3.1. Nature of the Host Galaxy 

The correlations between SFR and X-ray luminos- 
ity (Hornschemeier et al. 2005) and between Lhq and 
Lx(0.5 - 2keV) (Zezas 2000; Georgakakis et al. 2006) 
observed for star-forming galaxies would both imply 
Lx ~ 1 X 10"*° erg s^^ for D23H-1, which is consis- 
tent with the observed upper limit to the steady-state 
X-ray luminosity. If we were to assume that all of the 
Ha luminosity were instead powered by an AGN, which 
is not consistent with the diagnostic narrow-line ratios 
which indicate that the object is at most a compos- 
ite object, then the mean relation between Lho versus 
Lx{2 - lOkeV) for AGNs would predict a hard X-ray 
luminosity of L(2 — 10keV)« 10^^ erg s^^, which is a 
factor of 10 above the observed hard X-ray upper limit 
(Terashima et al. 2000; Ho et al. 2001). However, there 
is a significant scatter in the Lx versus Lho relation at 
lower luminosities {Lho < 10^^ erg s~^), and so we can- 
not rule out that the Ha line is being powered by an AGN 
from Lx/La alone. Finally, if the total [O HI] luminos- 
ity were powered by an AGN, then one would expect a 
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hard X-ray luminosity of > lO"'^ erg (Heckman ct 
al. 2005), which again is not observed. 

The lack of broad emission lines or high excitation 
narrow-line ratios in the host galaxy spectrum disfavor 
the presence of a normal AGN. However, UV variabil- 
ity has been associated with LINER galaxies (Maoz et 
al. 2005; Totani et al. 2005), and used as evidence that 
a low-luminosity AGN powers the low-ionization emis- 
sion lines, as opposed to a young stellar continuum. 
The largest amplitude NUV variability in these galax- 
ies reach factors of 3 — 5 (NGC 4552, Renzini et al. 1995; 
NGC 4579, Earth et al. 1996; NGC 1399, O' Connell 
et al. 2005; NGC 3368, NGC 3998, NGC 4203 Maoz et 
al. 2005), with luminosities of - lO^^^''^ erg s^^. Al- 
though the amplitude of NUV variability of these ob- 
jects are comparable to that observed for D23H-1, their 
UV fluxes were measured rising high spatial resolution 
Hubble Space Telescope (HST) imaging that isolates the 
nucleus from the host galaxy UV light. Our GALEX 
images in the NUV have a spatial resolution of 5", and 
can include UV starlight in the aperture. Given that 
we successfully fit the steady-state UV/optical SED with 
starlight alone, with no contribution from an AGN, it is 
likely that some if not all of the steady-state UV flux is 
from starlight. Thus, the observed variability amplitudes 
for D23H-1 in the FUV and NUV are merely lower limits 
to the intrinsic nuclear variability. 

BL Lacs a rare subclass of AGNs that have weak emis- 
sion lines, but are characterized by strong non-thermal 
continuum emission and rapid variability. Their pecu- 
liar properties are explained as the result of the orien- 
tation of a relativistic jet in these sources toward the 
direction of the observer (e.g., Blandford & Rees 1978). 
The difficulty with this scenario as an explanation for 
the lack of AGN spectral features in D23H-1, is that it 
is not detected in the FIRST radio catalog (White et 
al. 1997), with an upper limit to the flux at 1.4 GHz 
of 0.75 mjy. The corresponding radio-to-optical spectral 
index, a^o, defined by a^^^,^ = -log(/^i//i.J/log(i^i/i^2), 
where vi = 4400 A, and V2 = 1-4 GHz, is Q!ro=0.24, and 
places the galaxy in the radio quiet AGN category, and 
below the values seen even for BL Lacs with high-energy 
cutoffs (Bcckmann et al. 2003). Furthermore, the char- 
acteristic SED of a BL Lac has a parabolic shape that 
extends from the radio to the X-rays. The lack of a non- 
stellar continuum diluting the absorption features in the 
optical spectrum simultaneous with the peak of the flare 
is also inconsistent with this picture. Finally, the hosts 
of BL Lacs are typically massive elliptical galaxies, and 
not « 1O^"M0 star-forming galaxies. 

The UV/optical-to-X-ray spectral index observed dur- 
ing the flare, defined in the literature with vi = I keV or 
2 keV and 1^2 = 2500 Ais aox > 2.4 ± 0.2. Here, we are 
using the extinction-corrected NUV flux density f^^ = 
lQA^/2.5f^j(^l + z) ^ (5.6 ± 2.0) X 10-28 erg s^^ cm'^ 
Hz-i, where A(2316A/(1 + z)) = 3.95E{B - V)gas = 
1.2 ±0.4; and the extinction-corrected X-ray flux density 

hkeV{l+z) < 10^-/2-5/ife,v/(l + z) = (2.0 ±0.1) X lO--^-^ 

erg s-i cm-2 Hz-\ where 10^^/^-^ = (1.4 ±0.1). With- 
out the internal extinction correction of a factor of 3 ± 1 
in the NUV and 1.4±0.1 in the X-rays, Uox is > 2.2. The 
observed upper limit to the X-ray flux density is 2 orders 
of magnitude fainter than expected for Seyfert 1 galaxies 



in the same luminosity class, which have = 1-0 — 1.6 
(Strateva et al. 2005; Steffen et al. 2006) or BL Lacs 
{ctox = 0.6 — 1.6, Plotkin et al. 2008), but in the range ex- 
pected for normal galaxies {aox > !)• One can argue that 
the spectral index observed during a flaring state cannot 
be compared to the average spectral indices observed in 
AGNs. There may be a time delay between fluctuations 
in the UV band and the X-ray band, or spectral changes 
that occur in AGNs during high states of luminosity. It 
is a known property of AGNs that they have "bluer" col- 
ors when they brighten in the optical (Giveon et al. 1999; 
Geha et al. 2003). We also found the same behavior in 
the FUV— NUV colors during the high state of variable 
AGNs and quasars (Gezari et al. 2008). The implications 
of aox during the steady state are not clear, since we do 
not know what fraction (if any) of the steady-state NUV 
flux is from an AGN. If, for example, we assume that 
10% of the steady-state NUV flux is nonstellar, then aox 
of the nonstellar component is > 1.5, and is in the ob- 
served range for AGNs with !/L^(2500A)~ 10^^ erg s"^ 
(see Figure 2 in Maoz 2007). Thus, the UV and X-ray 
properties of the source during its steady state do not 
rule out the presence of a low-luminosity AGN. 

Considering that the host of the flare is a star-forming 
galaxy, we should also investigate the possibility that 
the transient source is a supernova. Hydrogen-rich core- 
collapse supernovae (Type II SNe) arc hot {Teff ^ 10^ 
K) and bright in the UV at early times (i.e. SN 1987A, 
Hamuy et al. 1988; SN 2006bp: Immler et al. 2007, 
Dessart et al. 2008; SN 2005cs, Dessart et al. 2008; SN 
2008es, Gezari et al. 2009). The SN ejecta cool through 
expansion and radiation, causing the peak of the SED 
to shift to optical wavelengths, and the UV emission to 
drop rapidly on the timescale of weeks. Although SNe 
Iln have a UV excess above the cooling blackbody emis- 
sion that is powered by interaction with circumstellar 
material that can last for years (i.e., SN 1993 J and SN 
1998S, Fransson et al. 2005), this UV emission is domi- 
nated by emission-line features in a nebular phase, and 
is a small fraction of the continuum emission from the 
ejecta. Thus, the timescale of the UV flare presented here 
of over a year, and the small change in the FUV— NUV 
color during that time, is not consistent with the rapid 
cooling associated with SN ejecta. These arguments in 
combination with the fact that the position of the source 
in the UV during all of the GALEX epochs is constant 
within the astrometric accuracy of 0'.'5 (for a galaxy at 
dL = 900 Mpc, 0'.'5 corresponds to 2 kpc), and that the 
UV source and the optical flaring source^ measured from 
the difference image are coincident with the centroid of 
the host galaxy favors the association of the flare with 
the central supermassive black hole. 

3.2. SED of the Flare 

We now analyze the SED of the flare, by subtracting off 
the "steady-state" FUV, NUV, and g band data points 

from the "flaring-state" flux, and fit the remaining SED 
with a single temperature blackbody in the rest frame 
of the galaxy. Assuming that the steady-state SED is 
dominated by star formation, we correct the fluxes for 
an internal extinction of E{B — V)gas = 0.3 using the 
Calzetti (2001) law. Since the flaring continuum is asso- 
ciated with the central black hole, the extinction in the 
line of sight to the galaxy nucleus may not be the same 
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as the global extinction to the galaxy measured from the 
Balmer decrement of the HII region emission lines. Thus, 
to be conservative, we examine both the extreme case of 
no extinction in the line of sight to the galaxy nucleus, 
together with the more probable case of an extinction 
equal to the global extinction in the galaxy. 

The least-squares fit to a blackbody assuming E{B — 
V)gas = 0.3 and the Calzetti et al. (2001) law, with 
internal extinction correction factors of 4.0, 3.0, 2.0 in 
the FUV, NUV, and g band, respectively, yields Tht = 
l.ltll X 10^ K and log(i6f,(erg s~i))=45.0j5^ The 
steady state SED fitted with the BC03 galaxy template, 
and the blackbody model fit to the flare added to the 
galaxy template are shown in Figure [7| in units of v^L^^ 
as a function of emitted frequency, v,, = z)vo, where 
— /y^47r(i|/(l + z). Alternatively, the power-law 
fit "to the UV/optical flare SED yields a = +1.8 ± 0.4, 
where ex i/", and L(1350 - 8500A)= 5.0 x 10''^ erg 
s~^. Also shown in Figure [7] is the case for no extinc- 
tion in the line of sight, which yields Thb = 4.91q;5 ^ 
K, log(L66(erg s-i))=43.4+°-2, or a = +1.2 ± 0.4 and 
L(1350 - 8500A)= 1.6 x 10''^ erg s"!. 

3.3. Nature of the UV/Optical Flare 

While there is no direct evidence for an active nucleus 
in the host star- forming galaxy of D23H-1, we cannot 
definitively rule out the presence of a low-luminosity 
AGN with Lx < lO'*^ erg s^^. There are two sce- 
narios associated with the host galaxy's central super- 
massive black hole (SMBH) that could produce the lu- 
minous (> lO''^ erg s~^) UV/optical flare: (1) a flar- 
ing low-luminosity AGN with peculiar X-ray and optical 
emission-line properties during its flaring state that differ 
significantly from measured average properties of AGNs 
and (2) the tidal disruption of a star by an otherwise qui- 
escent SMBH. The important distinction between these 
two scenarios is the presence or lack of a steady-state 
accretion disk around the central black hole. 

In the case of a tidal disruption event, the UV fiare 
is caused by the sudden increase in the supply of gas 
on to the central black hole from a stream of tidally 
disrupted stellar debris that forms a thick accretion 
disk or torus, and emits a soft, thermal spectrum with 
Te// « (Li5dd/47rcri?2,)i/4 (uimer 1999). In the ac- 
cretion disk UV flare is produced by an up- 
ward fluctuation in the accretion rate, possibly due to 
a magnetohydrodynamic (MHD) instability (see review 
by Schramkowski & Torkelsson 1996). The timescale of 
variability due to accretion disk instabilities has been in- 
vestigated (Kawaguchi et al. 1998; Hawkins 2002), and 
rest-frame NUV variability studies of the ensemble statis- 
tics of sparsely sampled optical light curves of tens of 
thousands of quasars are in general agreement with the 
slope of the structure functions in these models (Vanden 
Berk et al. 2004; Bauer et al. 2009). While a disk spec- 
trum during episodes of higher accretion rates is expected 
to have a higher effective temperature, and thus appear 
"bluer" (e.g., Shakura & Sunyaev 1973), there are few 
theoretical predictions for the corresponding changes in 
the X-ray properties of the disk. 

There are only a handful of AGNs for which simul- 
taneous UV and X-ray variability has been studied in 
detail (Uttley et al. 2003; Arevalo et al. 2008; Breedt 



et al. 2009), and they demonstrate a strong correlation 
between X-ray and optical flux variations on time-scales 
of days to months with no measured time-lag. Rever- 
beration mapping experiments of tens of nearby AGNs, 
meanwhile, measure a time lag between the UV/optical 
continuum and the fluxes of broad, recombination lines 
that reprocess that continuum on the order of weeks to 
months, corresponding to the light-crossing time of the 
broad-line region (Peterson et al. 2004). Future wide- 
fleld time domain X-ray surveys, such as EXIST (Grind- 
lay 2004), in coordination with the GALEX Time Do- 
main Survey in the UV, and optical synoptic surveys, 
such as Large Synoptic Survey Telescope (LSST), could 
place much needed observational constraints on the rela- 
tionship between optical, UV, and X-ray variability in a 
larger sample of AGNs. 

Given that both the UV/optical SED and the X-ray-to- 
UV flux ratio during the flare in D23H-1 are anomalous 
compared with the average properties of AGNs, we inves- 
tigate further the possibility that the flare was a result of 
the tidal disruption of a star by an otherwise quiescent 
central black hole, with the caveat that the temporal be- 
havior of these properties are not well constrained both 
theoretically or observationally in the alternative case of 
a large amplitude, high luminosity AGN flare. 

3.4. Light Curve 

The light curve of a tidal disruption event should fol- 
low the power-law decay of the mass accretion rate of the 
tidal disruption debris. The evolution of the mass accre- 
tion rate is predicted by analytical calculations (Phin- 
ney 1989) and numerical simulations (Evans & Kochanck 
1989; Ayal et al. 2000) to decline as [(t-ir,)/(to-iD)]~", 
where n = 5/3, to is the time of the peak of the flare, 
and t]j is the time of disruption. We perform a least- 
squares fit of a power-law decay to the FUV light curve 
of D23H-1, after subtraction of the steady state FUV 
flux. In the fit, to is constrained from the observations, 
to is unknown and is a free parameter in the fit, and the 
power-law index is fixed to n = 5/3. This fit yields a time 
of disruption of tjj = 2007.24+ 0.06. If we instead allow 
the power-law index to vary, we get to = 2007.4 + 0.4 
and n = 1.3 + 0.9. 

We can place a strict upper limit of to < 2007.67 from 
the UV detection before the peak of the flare, so we as- 
sume = 2007.4+°:f yr, which yields {to-tD)/{l+z) = 
0.06 — 0.6 yr. This time delay is related to the central 
black hole mass as Mbh = k^{Rp/ RT)~^r~^m1[{to — 
tD)/O.llyr]2M0 (Li et al. 2002), where k ranges from 1 
to 3 depending on how much the star is spun-up upon 
breakup, r^, = Ri,/RQ, and rrii, = M-^/Mq, and so 
we get a range of possible black hole mass solutions 
that depend of the spin-up and type of the star dis- 
rupted, and its penetration factor (Rp/Rx), A/bh = 
(0.3 - 30) x l(fMQk^{Rp/RT)-^r-^ml. 

In Figure [H we also show a flt to the analytical equa- 
tion from Lodato et al. (2009) (L09) for the mass ac- 
cretion rate of a tidally disrupted solar-type star with a 
polytropic index, 7 = 1.4, which takes into account the 
internal structure of the star in the calculation. The in- 
stantaneous power-law index of the calculated mass ac- 
cretion rate is < 1 at early-times, and asymptotically 
approaches n ~ 5/3 at late times. We take their equa- 
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tion for M{t) for solar-type star disrupted by a 10® Mq 
black hole at Rp — Rt, and scale it by A;~'^/^Afg^^(l + z), 
where Mg = AIbh/IO^Mq. There is a degeneracy in the 
fit between the mass of the black hole and the spin-up of 
the star, since decreasing the spin of the star or increas- 
ing the mass of the black hole both result in stretching 
out the light curve in time. If we assume that the star is 
spun-up to near break-up upon disruption, fc = 3, then 
Mbh = (5.4 ± 0.4) X 10^ Af©, and to = 2007.56 ± 0.01. 

3.5. Central Black Hole Mass 

The black hole mass inferred from the relation be- 
tween Mbh and Mbuige (Magorrian et al. 1998), using 
the galaxy mass derived from the galaxy template cor- 
responds to Mbh ~ 3 X 1O''(B/T)i-"M0, where {B/T) 
is the bulge-to-total luminosity ratio, which is unknown. 
For a star-forming galaxy with k, lO^^M© it is likely 
a spiral galaxy and {B/T) can range from 0.25 in Sb 
galaxies to 0.01 for Sd galaxies, or in the case of M33 
can be almost a pure disk galaxy {B/T < 0.01, van den 
Bergh 1991). Spiral galaxies also appear to have sys- 
tematically lower black hole masses for their given bulge 
mass compared with the Magorrian et al. (1998) relation 
(Salucci et al. 2000). Greene & Ho (2005) demonstrated 
that the velocity width of the [O IIIJA5007 emission line 
(cr[0 III]) can be used as a proxy for the velocity disper- 
sion of the stars in the host galaxy bulge (cr*). We find 
C[o III] — 140±20 km s^^ in the galaxy rest frame, where 

(T[o III] = V'^lbs ~ '^fnsf Using the Mbh - <y* relation 
from Merritt & Ferrarese (2001), and including the rms 
scatter in the log(cr[o m])— log((7i,) relation of 0.2 (Nel- 
son & Whittle 1996), this velocity dispersion corresponds 
to a central black hole mass of log(MBH/™sun) = (7±1). 
While the error bars of this estimate are large, it is en- 
couraging that the black hole mass range is within the 
masses for which the tidal disruption of a solar-type star 
is observable. 

Assuming that L{t) = io[(i — ^d)/(^o — ^_d)]~"^'^, where 
we set Lq to the lower-bound to the peak blackbody lu- 
minosity, Lq = 3 X lO^'' erg s~^, and io = 2007.74, and 
to = 2007.4, then the total energy radiated found by in- 
tegrating the curve from to out to infinity is 1 x 10^^ erg, 
and corresponds to a total material accreted of Mace = 
Etot/{^(?) = O.O6M0, for e = 0.1, and M* > O.IMq if at 
least half of the stellar debris is expelled after disruption 
and escapes the system. 

4. DEMOGRAPHICS OF TIDAL DISRUPTION EVENT 
CANDIDATES 

The value for the power-law index fitted to the light 
curve of D23H-1, n — 1.3 ± 0.9, has enough flexibility 
in its error bars to be within Ict of n = 5/3. The early 
light curves of the GALEX candidates Dl-9 and D3-13, 
detected in the g, r, and i bands by CFHTLS in the first 
4 — 5 months after the peak, were also well fitted with a 
n = 5/3 power law, however their late-time light curves 
traced by the UV photometry showed a break in the 
power-law to a shallower decay with njji^g = 1.1 ± 0.3 
and nD3-i3 = 0.82 ± 0.03 (Gezari et al. 2008). This 
behavior at late-times is actually a prediction of tidal 
disruption models, which predict that after the debris 
shocks and circularizes, the mass accretion proceeds 
through viscous processes which should have a charac- 



teristic power law of n = 1.2 (Cannizzo et al. 1990). The 
power law fits to the light curves of all three GALEX 
tidal disruption event candidates as well as fits to the 
analytically derived L09 curve for a solar-type star are 
shown in Figure El Assuming k = 3, the fits to the L09 
curve have Mbh = 2.4 x IO^A^q and A/bh = 1-2 x lO^Af© 
for Dl-9 and D3-13, respectively. To more clearly show 
the shape of the UV/optical light curves of Dl-9 and 
D3-13, the r-band data points are shown along with the 
NUV data, shifted up by their average NUV— r color: 
(NUV-r)£,i-9 = -2.0 and (NUV-r)D3_i3 = -1.4. The 
tidal disruption event candidates discovered in the X- 
rays by ROSAT and XMM Slew Survey (XMMSLl) have 
sparsely sampled X-ray light curves, but their dramatic 
fading 3 — 10 years after the peak of the flares are also 
consistent with an n — 5/3 power-law decay (Komossa 
et al. 2004; Halpern et al. 2004; Li et al. 2002; Esquej et 
al. 2008). 

It is surprising how well the light curves of the GALEX 
candidates are described by a simple power-law decay, 
and that they are consistent with a power-law index of 
n = 5/3 at early times. However, the value of the power- 
law index is very sensitive to the time of disruption, and 
vice versa. Thus if one wants to use the light curves to 
measure to — tjj and get an estimate the central black 
hole mass, or if one wants to measure the power-law in- 
dex, n, to understand the declining mass accretion rate 
as a function of time, it is more constraining to use mod- 
els that remove some degrees of freedom in the fit by 
predicting the rise of the light curve as well as its decay, 
such as those calculated by L09. 

Shown in Figure [9] are the blackbody fits and the 
power-law fits to the UV/optical SEDs of the three 
GALEX candidates, and in Table [2] we list the fitted pa- 
rameters, the reduced statistic {xu), and the resulting 
luminosities. Because the peak of the flares from Dl-9 
and D3-13 were detected in the optical, with the first 
UV detection ~four months after the peak, we estimate 
the simultaneous UV/optical SED by extrapolating the 
UV fluxes to the time of the peak using the power-law 
flt to the light curves (shown in Figure [8]). We also show 
the single blackbody flts to the soft X-ray detections for 
Dl-9 and D3-13, that were detected 2.3 and 1.2 yr af- 
ter the peak, respectively, with Tbi, ~ 2.7 x 10^ K and 
Tfcb = 4.9 X 10^ K, respectively, which yield blackbody 
luminosities of 5.7 x lO''^ erg s~^ and 2.3 x 10''^ erg s~^. 
It is important to note that these X-ray detections did 
not have enough counts to measure the spectrum or its 
associated blackbody temperature directly, however, the 
lack of counts > 1 keV in the sources do place constraints 
on a power-law spectrum, fx , to have F > 3 and 
F = 7 ± 2, respectively, shown in Figure [D 

The power-law flts to the UV/optical SEDs of all 
three GALEX candidates are steep (a = -1-1.4 to +1.8), 
and not typical of an AGN power-law spectrum with 
-2.0 < a < -0.43 (Zheng et al. 1997; Dietrich et 
al. 2002). A single blackbody model appears to pro- 
vide good fit to the UV/optical part of the fiare SED 
of Dl-9 (x^ = 0.97) and for the internal extinction cor- 
rected fiare SED of D23H-1 {xl = 0.98). The poor flt to 
the UV/optical SED of D3-13 with a single temperature 
blackbody would be greatly improved with the addition 
of a second lower temperature blackbody to account for 
the flattening of the optical SED slope. However, the 
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blackbody model fits to the UV/optical SEDs need not 
be unique. One could imagine that the UV emission 
is the result of inverse Compton scattering from a hot 
plasma (corona) with seed photons from synchrotron ra- 
diation (perhaps in the radio or infrared). Future radio 
and infrared observations would be critical for testing 
this alternative scenario. 

Now that the census of tidal disruption candidates has 
reached a total of 10, we can start to compare their prop- 
erties and investigate the selection effects that arise from 
their different modes of discovery. The RASS detected 
large amplitude, luminous soft X-ray outbursts from four 
galaxies with no previous evidence of AGN activity. The 
soft X-ray flares were well described by blackbody spec- 
tra with Tbb = (6 — 12) X 10^ K, although they could also 
be fitted with extremely soft power laws with F = 3 — 5 
(Bade et al. 1996; Komossa & Bade 1999; Komossa & 
Greiner 1999; Grupe et al. 2003; Greiner et al. 2000). 
More recently a flaring soft X-ray source was reported 
from a member of the cluster Abell 3571 with a peak 
luminosity of L(0.3 - 2.4 keV) = 6.8 x 10'*^ erg s'^ 
and blackbody temperature of kTtb = 0.12 keV, which 
decayed by a factor of >650 over 10 yr (Cappelluti et 
al. 2009). The XMM Slew Survey (XMMSLl) detected 
large amplitude soft X-ray flares from two galaxies with 
no previously known AGN. Only one of the flares had 
enough counts to be directly fitted with a blackbody with 
Tf,b = 1.1 X 10^ K or a power law with F = 3, while the 
other had counts only in the 0.2 — 2 keV band, and was 
assumed to have T;,f, = 8 x 10^ K (Esquej et al. 2008). 

All of the candidates mentioned above have quiescent 
galaxy hosts, with the exception of NGC 3599 and NGC 
5905. NGC 5905 is a ROSAT candidate that is a star- 
burst galaxy that was found with narrow-slit STIS spec- 
troscopy to host a low-luminosity Seyfert 2 nucleus with 
Lhc 10^^ erg s-i (Gezari et al. 2003). NGC 3599 is 
an XMMSLl candidate for which follow-up spectroscopy 
revealed narrow emission lines with ratios on the border- 
line between a LINER and a Seyfert galaxy (Esquej et 
al. 2008). The presence of a low-luminosity AGN in these 
galaxies does not contradict the tidal disruption scenario, 
since their amplitudes of ^100 in the soft X-rays are 
more extreme than any previously known cases of AGN 
variability. 

Similarly to NGC 5905, D23H-1 could also harbor a 
Seyfert 2 nucleus that is masked by the surrounding HII 
regions in its spectrum. However, the hard X-ray lumi- 
nosity upper limits during all three GALEX flares place 
an upper limit on the presence of a persistent power-law 
AGN down to L{2 - 10) keV < 10^^ erg s'^ in D3-13 
and D23H-1, and down to L(2 - 10keV)< 6 x lO''^ erg 
in Dl-9. Although these X-ray luminosity upper lim- 
its do not exclude the presence of a low-luminosity AGN 
in the nuclei of the galaxies, the unique SEDs and light 
curves of the flares make them stand out from known 
low-luminosity AGN behavior. 

If we interpret the steep power-law shapes of the UV 
and soft X-ray candidates as the Rayleigh Jean's and 
Wien's tail of a thermal blackbody spectrum, respec- 
tively, then their blackbody temperatures, in combina- 
tion with their bolometric luminosities can be used to 
infer a characteristic radius of the emission. Figure [10] 
shows the peak bolometric luminosity and blackbody 
spectrum temperature {kT^b) of the four candidates from 



RASS summarized in Komossa et al. (2002), the Abell 
cluster candidate from Cappelluti et al. (2009), two can- 
didates from XMMSLl reported by Esquej et al. (2006, 
2008), the candidate presented in this paper, and the 
two previous candidates discovered by GALEX reported 
by Gezari et al. (2006, 2008). The blackbody compo- 
nent of D23H-1 is shown both without the internal ex- 
tinction correction (solid symbol) and with a correction 
for E{B — V)gas = 0.3 (open symbol), and connected 
with a solid line. For the GALEX candidates that had 
a soft X-ray detection during their flare, we connect 
the UV/optical blackbody component to the soft X-ray 
blackbody component (plotted with an X) with a dotted 
line. Because the X-ray detections for these candidates 
were 1 — 2 yr after the peaks of the flares, we scale 
the fluxes using the light curve fits to estimate their X- 
ray luminosities during the peak. However, this is highly 
uncertain, since the X-ray detection for D3-13 showed 
variability of a factor of 10 on a 1 day timescale, and 
may not follow the decay of the UV/optical light curve. 
We plot the soft X-ray component of D3-13 as an upper 
limit to account for this variability. 

Also shown is the UV flare detected by iJST/Faint Ob- 
ject Camera (FOC) in the nucleus of NGC 4552 (Renzini 
et al. 1995; Cappellari et al. 1999) interpreted as the ac- 
cretion of a tidally stripped stellar atmosphere by the 
mini- AGN at its center. This galaxy was shown to have 
another episode of variability 10 yr later, when its NUV 
flux increased by a factor of 1.2 between epochs separated 
by 2.4 months (Maoz et al. 2005). This lower-amplitude 
variability most likely originates from the low-luminosity 
AGN detected in the nucleus of the galaxy via broad op- 
tical emission lines and a hard X-ray point source with 
4 X 10^^ erg s"^ (CappeUari et al. 1999; Xu et 
al. 2005). Given the rarity of tidal disruption events, 
the fact that several other nearby LINERs demonstrated 
UV variability of a similar amplitude as the flare in NGC 
4552 reported a decade earlier by Renzini et al. (1995) 
casts doubt on its interpretation as a tidal stripping 
event. 

Shown in comparison are the luminosities expected for 
a blackbody radius equal to the tidal disruption radius 
of the central supermassive black hole, using Lboi = 
4:TrRl^a{Tbb/ fc)'^, where fc = Tbh/T^/f takes into ac- 
count deviations of the color temperature from the ef- 
fective temperature due to electron scattering and comp- 
tonization in systems with high-mass accretion rates (see 
Li et al. 2002). Also shown is the Eddington lumi- 
nosity of black holes of 10^, 10^, and IO^Mq, where 
^Edd = 1-3 X 10^^ Mg erg s~^. There is a large range 
in color temperatures of the events, with a strong bias 
toward the wavelength range of the survey used to dis- 
cover them. The temperatures do not correlate with the 
luminosity, as would be expected if there were a single 
characteristic radius from which the emission was pro- 
duced. 

The black hole masses inferred from the host galax- 
ies of the candidates using locally established scaling re- 
lations between A/bh and the host galaxy bulge, yield 
black hole masses for the XMMLSSl candidates of a 
few times lO'^M©, ~ lO'^M© for the AbeU cluster can- 
didate (Cappelluti et al. 2009), and a few times lO^M© 
for the GALEX candidates and the ROSAT candidate 
NGC 5905. The low luminosity (- lO^^LEdd) of NGC 
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5905 was explained by Li et al. (2002) as a result of 
the disruption of a brown dwarf or giant planet, while 
they attributed the small inferred radius as emission 
from a small patch possibly originating from a nar- 
row stream of accreting material, instead of an opti- 
cally thick circular torus of material at the tidal dis- 
ruption radius with T^fj « (LBrfrf/47rcri?|,)^/* = 2.5 x 

W^KM^^^'^r'^^^m^^^^ (Ulmer 1999). 

In contrast, the blackbody radii inferred for the 
UV/optical component of the GALEX flares are 10 times 
the tidal disruption radius of a IO^Mq black hole, shown 
with a dotted line in Figure [TUl The presence of material 
at such large radii is a natural consequence of the dis- 
ruption of a star, since at least half of the stellar mass is 
unbound at high velocities during the disruption process. 
An extended envelope of debris may form and reprocess 
the radiation from the inner accretion torus with an effec- 
tive temperature of Te// « 2.3 x lO^K {Mbh/W'^ Mq)^/'^ 
(Loeb & Ulmer 1997; Ulmer et al. 1998), close to the 
temperatures fitted to the UV/optical components of the 
GALEX flares. 

The events, regardless of what wavelength range they 
were discovered in, share the properties that (1) their 
peak accretion rates are large enough to produce high 
luminosities (10^^ erg s""'^), (2) their SEDs appear to be 
steep power-laws consistent with being thermal in na- 
ture, and (3) the power-law decay of their light curves 
is consistent with the declining mass accretion rate ex- 
pected from the fallback of debris from a tidally disrupted 
solar-type star. However, the geometry of their radiat- 
ing regions is diverse, and there appear to be two com- 
ponents to their continuum emission. The soft X-ray 
components have radii on the order or smaller than the 
tidal disruption radius of the central black hole, and are 
most likely associated with the inner accretion torus (or 
a small region of it) . The rapid variability of the soft X- 
ray source detected in D3-13 (Gezari et al. 2006, 2008) 
is consistent with this picture. While the UV/optical 
components of the flares appear to have a characteristic 
radius of roughly 10i?T, possibly tracing an envelope of 
debris surrounding the black hole. 

5. CONCLUSIONS 

We present the third candidate tidal disruption event 
discovered by GALEX: a UV/optical flare from a star- 
forming galaxy aX z — 0.1855 with a broadband SED, 
light curve, and peak luminosity in excellent agreement 
with the expected properties of a flare from the tidal dis- 
ruption of a star. While we can only place upper limits on 
the presence of a steady-state low-luminosity AGN in the 
host galaxy down to Lx < 10**^ erg s"-^, the broadband 
properties of its flare deviate from the average properties 
observed for AGNs. While it is possible that the flares we 
are detecting with GALEX represent a large amplitude, 
high-luminosity tail of the distribution of UV variability 
from AGNs, such flares would have to comprise of UV 
emission that is not correlated with hard X-ray emission 
or optical emission lines (unlike what is observed in the 
small number of AGNs that do have detailed optical and 
X-ray monitoring campaigns). Future coordinated opti- 
cal, UV, and X-ray monitoring of AGNs over wide flelds 
of view will be necessary to determine how favorable this 
interpretation is over the tidal disruption event scenario. 

The next generation of wide-fleld synoptic surveys will 



be at optical wavelengths, and their capabilities for de- 
tecting tidal disruption events will depend on how these 
flares manifest themselves in this wavelength range. We 
have now successfully detected three candidates for tidal 
disruption events from quiescent SMBHs in the optical, 
and in Gezari et al. (2008) we made favorable predictions 
for the detection rates for Pan-STARRS and LSST us- 
ing the theoretically derived black hole mass-dependent 
volume rate of tidal disruption events that successfully 
reproduced the observed detection rate by GALEX in 
the UV. For this calculation, we assumed that flares at 
their peak radiated at the Eddington luminosity with 

Tbb — 2.5 X lO^Afg^^ K, but that the fraction of events 
that produced Eddington luminosity flares declined with 
increasing black hole mass (Ulmer 1999), and we assumed 
the light curves followed power-law decays with n = 5/3, 

^0 ^ = 0.11 yr fc^'^/^Mg^^, and a spin-up parameter 
for the star oi k = 3. 

Here, we will use a simpler approach to estimate 
the detection rate of tidal disruption events in opti- 
cal surveys by using the observed range of Galactic 
extinction-corrected peak absolute g-band magnitudes 
for the GALEX candidates, Mg(Dl-9) = -17.66±0.05, 
Mg{D3 - 13) = -18.88 ± 0.05, Mg{D23H - 1) = 
-17.52 ± 0.12, where Mg at peak for Dl-9 and D3- 
13 is estimated from the peak observed r-band magni- 
tude and the average observed g — r color during the 
flares. We estimate the black hole mass-dependent vol- 
ume rate as iV(MBH) = A^(MBH)A^Bff (Mbh)/(Mbh), 
where N{Mbh) is the tidal disruption rate as a function 
of black hole mass found for elliptical galaxies by Wang 
& Merritt (2004), 7V(Mbh) ~ 1-6 x lO"-* yr-lMg-°■^ 
Nbh is the black hole mass function estimated from the 
Ferguson & Sandage (1991) E-I-SO luminosity function 
scaled to black hole mass using the A/bh — -Aibuigc re- 
lation from Merritt & Ferrarese (2001) and the mean 
galaxy mass-to-light ratio from Magorrian et al. (1998) 
with a factor of 2 to account for the bulges of spiral 
galaxies, and /(Mbh) = -0.9091og(MBH/MQ) + 7.00 is 
the fraction of flares that have a mass accretion rate close 
to the Eddington limit for a Salpeter mass function from 
Ulmer (1999). The declining fraction of Eddington lumi- 
nosity flares from higher mass black holes is a result of 
the decrease in peak accretion rate with increasing black 

— 1/2 

hole mass as . If we integrate TV (Mbh) over the 

full range of black hole masses, 10^ — IO^A/q, this yields 
iV = 1.3 X 10-5 yr-i Mpc'^. 

For this volume rate, the 50 deg^ Pan-STARRSl 
Medium Deep Survey (PSl MDS) will detect 32 events 
yr-i with 5 < 22 out to 790 Mpc (or z = 0.16). We 
have chosen this brightness cut off so that the peak of 
the flares is at least 2 mag brighter than the detection 
limit of the survey to allow for the decay of the light 
curve of the flare to be followed in detail. This detection 
rate is double our previous estimate of 15 events yr"^. 
However, here we have assumed that all flares will have a 
peak with Mg < —17.5. If instead we take into account 
that these luminosities apply only to the higher mass 
black holes, i.e., AfsH > lO^Af©, then the volume rate 
for this black hole mass range, 2.3 x 10~^ yr""'^ Mpc~'^, 
corresponds to a detection rate for PSl MDS of 5 yr~^. 
This can be thought of as a lower limit, since the survey 
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will be sensitive to fainter events from lower mass black 
holes at smaller distances. The future LSST will act as 
a tidal disruption event factory, since its depth and ca- 
dence will be similar to PSl MDS, but with a factor of 
400 increase in area that will increase the detection rates 
by 2 orders of magnitude! 

Shallower, wide-field optical time domain surveys such 
as the PSl Stt survey, the Palomar Transient Factory 
(Rahmer et al. 2008) and Skymapper (Keller et al. 2008) 
have the potential to detect bright events in the local uni- 
verse, with detection rates per year of w 1 event (1000 
deg^)~^ with g < 19 out to 200 Mpc, or w 1 event 
(1000 deg^)~^ with r < 19, assuming g — r > —0.4, 
as observed for the GALEX candidates Dl-9 and D3-13, 
{g - r)Di-9 = -0.1 and {g - r)D3-i3 = -0.4. 

In the near future, we will have the opportunity to ex- 
ploit the sensitivity, large volume, and daily cadence of 
the upcoming optical synoptic surveys to accumulate a 
large number of detailed light curves of tidal disruption 
events in the local and distant universe. While each in- 
dividual tidal disruption event provides a cosmic labora- 
tory to the study the detailed accretion physics around 
supermassive black holes, the ensemble properties and 
rates of a large number of events will bo a powerful tool 
for testing the predictions of galaxy stellar dynamical 
models, and probing the mass function of dormant su- 
permassive black holes in normal galaxies. 
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Fig. 1. — GALEX yearly epochs of observations that show the large amplitude flare with a factor of 9.0 increase in flux in the FUV and 
factor of 4.7 increase in flux in the NUV between 2006 and 2007. Circles show the circular aperture with r = 6'.'0 in the FUV and r = 3'.'8 
in the NUV. 



TABLE 1 
GALEX Observations 



UT Date texp (ks) FUV NUV 



2004.798 


17.3 


22.81 


± 


0.13 


22.04 


± 


0.20 


2005.753 


6.3 


22.93 


± 


0.21 


21.98 


± 


0.20 


2006.693 


12.3 


23.10 


± 


0.16 


22.30 


± 


0.20 


2007.673 


16.1 


21.24 




0.07 


21.18 


± 


0.19 


2007.744 


6.1 


20.71 


± 


0.09 


20.61 


± 


0.19 


2007.794 


7.4 


20.86 


± 


0.06 


20.79 


± 


0.20 


2008.685 


10.3 


22.30 


± 


0.09 


21.86 


± 


0.20 


2008.724 


10.6 


22.14 


± 


0.10 


21.82 


± 


0.20 


2008.830 


17.6 


22.18 


± 


0.08 


21.82 


± 


0.20 
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2009 



Fig. 2. — Left: FUV and NUV GALEX light curve, and FUV-NUV color as a function of time, in AB magnitudes corrected for Galactic 
extinction. The source is bluer during the "flaring state" , indicating that the flare has a different SED from the "steady state" . 




Fig. 3. — Optical imaging of the field of D23H-1. Left: SDSS 3-band reference image, taken on 2001 October 15 UT. D23H-1 is indicated 
with the black tick marks. Center: LFC (/-band image of the same field on the night of 2007 October 15. Right: resulting subtraction, 
showing a residual point source at the location of the centroid of D23H-1. The two poorly subtracted stars were saturated in the LFC 
image. All images are oriented with north to the right and east up, as indicated by the compass. 
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Fig. 4. — Spectra from LRIS taken 13 days before the peak of the flare, and with DEIMOS 8 days after the peak of the flare. The 
best-fitting Bruzual & Chariot (2003) galaxy template is overplotted in green. Emission-line features are labeled with the solid tick marks, 
and stellar absorption features are labeled with the dotted tick marks. The hatch marks indicate the locations of the strong O2 telluric B- 
and A-band absorption features. 



TABLE 2 
UV/Optical SED Fits 



ID 


E(B - V),„t 


Tbb 


Tbb - lo- 


Tbb + lo- 


xUbb) 


a 


xUpi) 


log(i6b) 


log(Li,6 - lo-) 


log(L6b + lo-) 


log(ipi) 






(xlO-* K) 


(xlO* K) 


(xlO* K) 








(ergs s~i) 


(ergs s"-"^) 


(ergs s"-"-) 


(1350 - 8500 


Dl-9 


0.0 


5.6 


4.4 


6.1 


0.97 


1.4 ± 0.3 


1.77 


43.6 


43.5 


43.7 


43.3 


D3-13 


0.0 


4.9 


4.7 


5.1 


22.48 


1.1 ± 0.1 


14.50 


43.8 


43.7 


43.9 


43.5 


D23H-1 


0.0 


4.9 


4.4 


5.6 


0.64 


1.2 ± 0.4 


3.47 


43.4 


43.3 


43.6 


43.2 


D23H-1 


0.3 


17 


12 


39 


0.98 


1.8 ± 0.4 


3.12 


45.5 


45.0 


46.8 


43.7 
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Fig. 5. — Gaussian fits to the narrow emission lines (shown in blue) after the galaxy template (shown in green) has been subtracted 
account for the contribution of stellar absorption lines. 
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Fig. 6. — BPT diagram of diagnostic line ratios which classify D23H-1 as a star-forming galaxy. The empirical dividing line between 
star-forming galaxies and AGNs from Kauffmann et al. (2003) shown with a dashed line. Definitions of HII galaxies, LINERs, and Seyfert 
galaxies from Ho et al. (1997) shown with the solid lines. The dotted lines show the theoretical divisions from Kewley et al. (2001). 
The region between the dotted lines in the first panel is where composite objects lie. The gray dots show the "mixing line" between a 
pure star-forming spectrum and a LINER spectrum, by adding a contribution from LINER emission to the Balmer emission lines in 10% 
increments. 
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Fig. 7. — SED of the source in the "steady state" (black symbols) and in the "flaring state" (blue symbols) as a function of frequency 
shown in the rest-frame of the host galaxy at 2 = 0.1855. The fluxes have been corrected for Galactic extinction of E{B — V) = 0.035 
and extinction internal to the host galaxy of E{B — V)gas = 0.3 using the Calzetti (2001) law. The best-fitting Bruzual &; Chariot (2003) 
galaxy template to the "steady-state" SED is shown in red. A blackbody spectrum has been added to the host galaxy template to match 
the "flaring-state" FUV, NUV, and g-band flux with Ti,i, = 1.7 X 10^ K (thin line) when flare fluxes are corrected for internal extinction 
(open symbols), or Ti,i, = 4.9 X 10** K (thick line) when not corrected for internal extinction (filled symbols). The Chandra upper limits on 
the X-ray fiux densities at 1 keV are plotted with an X. 
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Fig. 8. — Left: Light curves of the three GALEX tidal disruption flare candidates, presented in this Paper and Gezari et al. (2008). The 
AB magnitudes have been corrected for Galactic extinction, and subtracted by the distance modulus (DM=5 log(dj^/10pc)). The time axis 
is in years since the time of disruption found in the fit, corrected for time dilation, in the rest frame of the galaxy. The FUV light curve is 
shown for D23H-1, and the NUV light curves are shown for Dl-9 and D3-13. The r-band light curves are shifted up by the observed average 
(NUV-r) color, (NVY-r) = —2.0 and {NUV-r)j33_i3 = —1.4. The solid lines show the power-law fit with the power-law index, n, 
flxed to 5/3, and the dashed lines show the fits with the power- law index n free to vary. Right: the solid line shows the analytical solution 
from Lodato et al. (2009) of M{t) for a solar-type star with a 7 = 1.4 polytropic index, scaled in time by k~^^'^(M-gii/10^ Mq)^^"^ , where 
we assume fc = 3, and fit for A/bh = 5.4 X 10^^,2.4 X 10^, and 1.2 X lO'^M© for D23H-1, Dl-9, and D3-13, respectively 




Fig. 9. — Left: single temperature blackbody fits to the three GALEX fiares. Fits to the soft X-ray fiux densities in Dl-9 and D3-13 are 
shown with a dashed line, and were detected 2.3 and 1.2 yr after the peak of their flares, respectively. The fluxes corrected for an internal 
extinction of E{B — V)gas = 0.3 using the Calzetti (2001) extinction law are plotted with open symbols, and the corresponding blackbody 
fit is plotted with a dotted line. Right: power-law fits to the three GALEX flares and the soft X-ray detections in Dl-9 and D3-13. The 
fluxes corrected for an internal extinction are plotted with open symbols, and the corresponding power-law fit is plotted with a dotted line. 
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Fig. 10. — Peak bolomctric luminosity of candidate tidal disruption events as a function of their blackbody temperature. GALEX 
candidates with soft X-ray detections during their flares have the estimated peak soft X-ray blackbody component plotted with an X, and 
connected to the UV/optical blackbody component with a dotted line. The internal extinction corrected blackbody component of D2.3H-1 
is plotted with an open symbol, and connected to the uncorrect blackbody component with a solid line. The solid gray lines show the 
luminosity as a function of temperature for a blackbody radius equal to the tidal disruption radius of a 10^, 10^, and 10** Mq black hole, 
and the dotted gray line shows the luminosity for blackbody radius equal to 10 times the radius of a IO^A/q black hole. The dashed gray 
lines show the same relation, but for T^b/T^ff = 3. The dashed black lines show the Eddington luminosity of a 10®, lO'^, and 10* Mq black 
hole. 



